function R = MakeLab8Report(Dfx, a, b, y0, eps, n0)
R = cell(3, 1);
[P, n, t] = EulerMethodWithInitialStep(Dfx, y0, a, b, n0, eps);
PlotPartialFunction(P);
h = (b - a) / n;
R(1) = cellstr(sprintf('Euler method: h = %.4f ; time = %.4f', h, t));
[P, n, t] = RungeKuttaMethodWithInitialStep(Dfx, y0, a, b, n0, eps);
h = (b - a) / n;
R(2) = cellstr(sprintf('Runge-Kutta 2nd order method: h = %.4f ; time = %.4f', h, t));
[P, n, t] = AdamsMethodWithInitialStep(Dfx, y0, a, b, n0, eps);
h = (b - a) / n;
R(3) = cellstr(sprintf('Adams method: h = %.4f ; time = %.4f', h, t));
end

